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An effective mathematical framework based on Presymplectic Geometry for dealing with the 
"phase space picture" of timeless dynamics in General Relativity is presented. In General Relativity, 
the presence of the scalar Hamiltonian constraint which vanishes leads to the problem of time, which 
can be solved, up to an extent by adopting a timeless formalism. This has been done by Carlo Rovelli 
and Julian Barbour et al. In this paper we present a phase space reformulation of Barbour's theory. 
' The Presymplectic dynamics of general totally constrained, reparametrization invariant theories is 

y—( , developed, then applied to Jacobi mechanics, relational particle mechanics and the dynamics of the 

f"^ ■ free relativistic particle. 
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I. INTRODUCTION 



o 

The problem of time has been one of the most extensively investigated topics in quantum gravity. Due to its key 
| | role, its implications have been studied from philosophical, conceptual, structural and mathematical points of view. 
bsO Not only has the foundational aspects of canditate quantum theories of gravity been much clarified but also the 
possible meaning of time itself, especially at a fundamental level. Of the several approaches that have come to the 
. forefront of curent research, two are markedly similar. One, pioneered by Rovelli, seeks to "forget time" conceptually. 
The other, initiated and developed by Barbour, does away with time from the very outset. This latter approach has 
recently been the topic of extensive investigation by a host of researchers. However, a unified mathematical formalism 

■ that embodies these ideas would be instructive and desirable. The presymplectic framework presented in this paper 
seeks to do exactly this. In the first of a two-part paper we develop this formalism in the context of general totally 

■ constrained, reparametrization invariant theories. We then apply it to Jacobi mechanics, relational particle me- 
0^ ' chanics and the free relativistic particle. The application of this to the case of gravity is treated in Part 2 of this paper. 

o: 

This Paper is organized as follows. In Section II we very briefly state the problem of time. In Section III, we 
discuss the Jacobi-Barbour-Bertotti theory. In Section IV we make some remarks on Ephemeris time and its role. In 
J> ■ Section V we take a look at Rovelli's theory. In Section VI we discuss the violation of Dirac's theorem for completely 
constrained systems. In Section VII we treat phase space as a presymplectic manifold. Section VIII presents the 
presymplectic formalism. Section IX carries some concluding remarks. 



II. THE PROBLEM OF TIME IN QUANTUM GRAVITY 

The problem of time in quantum gravity is perhaps as old as the subject itself. It was noted first by Bryce de 

Witt in his paper " The Quantum Theory Of Gravity" |4j in which he presents the constraint equations for canonical 
Quantum General Relativity, 

7T* = (2.1) 

= (2.2) 

"H* = (2.3) 

X 4 * = (2.4) 
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where the first two constraints correspond to the vanishing of momenta associated to the Lapse and the Shift, while the 
other two are the scalar Hamiltonian constraint and the vector Diffconiorphism constraint respectively. He then notes 
that in such a setting, if we were to consider the wave function <J> as a function ^(x, x°) of the spatial coordinates and 
time, we see that the expectation value of the field operator, let's say, the momentum conjugate, undergoes something 
peculiar: <fV J (x, 0)* = (x,x )^ 
This leads him to conclude 

" Since the statistical results of any set of observations are ultimately expressible in terms of expectation 

values, one therefore comes to the conclusion that nothing ever takes place in quantum gravity dynamics, 
that the quantum theory can never yield anything but a static picture of the world." 

Thus the Wheeler dc Witt equation, in the standard interpretation leads to the " frozen formalism" which apparently 
tells that nothing should happen in Canonical Quantum Gravity. So how can one attain at least a partial solution to 
this problem? In order to answer this, we must first note that the problem arises fundamentally due to the fact that 
the Hamiltonian is but a sum of constraints, and that they in themselves saturate the theory. Due to the absence of 
the unconstrained Hamiltonian, we are left with no operator which can generate infinitesimal time translations of the 
system. 

Now here's the question we shall attempt to answer in the following part of this paper and the next: How can we 
extract dynamics from a system that is totally constrained? We shall adopt the "No Time" interpretation of Barbour 
et. al. for both the present and the forthcoming models including canonical general relativity. 



III. JACOBI-BARBOUR-BERTOTTI THEORY 



We shall first attempt to understand what it means to eliminate "time" from the action. The essence of doing so 
would be to replace time by an unphysical evolution parameter, with respect to which we can differentiate coordinates 
to attain velocities, which can indeed exist with the notion of an absolute or even a clock time, and the action should 
really be invariant under an arbitrary change of this parameter, i.e be reparameterization invariant. Thus eliminating 
time really amounts to getting rid of a preferred time scale, i.e. the system should refer to no external clock. As we 
shall see in the coming chapters, in doing so, the system finds a way to be its own clock! Now if we considered A to 
be the evolution parameter, this condition implies that A — > /(A) should leave the action invariant. As we shall see, 
the Action for an ordinary nonrelativistic particle dynamics just wont do. This is because, 

S' = Jd\f'(\)(j±^T-U) (3.1) 

where, all the primes and velocities here denote derivatives with respect to A). 

As it is clearly evident, the Action is not invariant under such a transformation, largely due to the dependence of 
the Lagrangian on the square of the first time derivative of the position. So in order to attain a reparameterization 
invariant Action, we need to find an Action that is linear in velocities, so that we can perform a Routhian reduction 
by letting the Newtonian time 't' be a function of A and eliminating it. In doing so we obtain the Action 

Sjbb = J d\2y/(T(E-U)) (3.2) 

which is the Jacobi Barbour Berotti action that obeys the Jacobi geodesic variational principle, which tells us that 
the geodesies on configuration space extremize the Action. 



IV. EPHEMERIS TIME 



As mentioned earlier, for the system we just considered, the removal of time scale via reparameterization invariance 
of the Action would lead to an interesting result. The system becomes its own relational "clock", i.e. the notion of 
time is derived from dynamics. In order to exhibit this, we shall use a rather pedestrian derivation of the "Ephemeris 
Time" . We consider 

T = E-U (4.1) 
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where, the rj is the Euclidean metric. 

What we have here, in itself is a remarkable result, it is as though an increment of time is the "mass weighted 
hypotenuse" of the total dynamical increment in the system. Now if we make our earlier assumption that velocities 
should be defined over an arbitrary evolution parameter, we see that the above expression becomes, 

dt = {^jdX (4.4) 

=>h-t x = [{^jdX (4.5) 

Thus it is evident that the Ephemeris gives us a metric-trans-temporal notion of identity between two subsequent 
configurations, and time is now completely emergent from the dynamics. We shall not take a look at Rovelli's theory. 



V. ROVELLI'S THEORY 



Here we shall breifly discuss Carlo Rovelli's timeless theory. We shal be using the same mathematical framework 
developed by Rovelli et. al. but our formalism differs subtly in interpretation and mostly in application, as we use 
this formalism via the Barbour interpretation and not the evolving constants of motion. In this theory, the Phase 
space T*C 3 (q a ,P a ) is a 2n dimensional manifold that admits a k dimensional surface £ and carries the canonical 
one form 

0= Pa dq a (5.1) 

and 9\z =8, now, the dynamics is governed by the equation 

d6{X) = (5.2) 

where, X belongs to the 2n — k dimensional kernel of d9. Now the motions are integral surfaces of this equation. The 
Hamilton- Jacobi formalism of this theory is as follows. The Hamilton- Jacobi equation is given by 

where, S{q a ) is Hamilton's principal function. This admits a one parameter family of solutions S(q a , Q l ) parameterized 
by n — k constants of motion Q l and the equation for the "motions" of the system is given by, 

f( q «,Q\p i ) = — -P i =Q (5.4) 
for n — k constants Pi. Q, P coodinatize the 2(n — A:)dimensional phase space. 



VI. VIOLATION OF DIRAC'S THEOREM 



It is instructive to note here the invalidity of Dirac's theorem regarding first class constraints for completely con- 
strained Hamiltonian systems. 6Sy stems with first class constraints We begin by considering constrained Hamiltonian 
systems with first class constraints. The Hamiltonian is given by 

H T = H + v a <t> a (6.1) 

where the v a are arbitrary Lagrange multipliers and the </> are the primary first class constraint functions. For some 
function / which is a dynamical variable, an infinitesimal increment in time 8t is given by 

Sf = 5t{f,H T } (6.2) 
= 5t({f,H} + v a {f,<j> a }) (6.3) 

Due to the arbitrariness in the v a , we can alternatively choose v' a and the difference 



Af(St) = e a {f,cp a } 



(6.4) 
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where. 

e a = St(v a - v' a ) (6.5) 

is very small but due to the fact that such changes of v a should not lead to physically inequivalent motions, it is 
concluded that cf) a generates unphysical gauge tansformations. We must note that Dirac assumes that St is the same 
for both v a and v' a . This is the crux of the matter as it tells us how absolute time influences this interpretation. 
BRepametrization invariant models Here the Hamiltonian is a constraint 

H T = V a (f>a (6.6) 
an example of which is the Barbour-Bertotti theory where 

H ^N{^-+V -E) (6.7) 

777 

We note that in this model there is no preferred notion of time. Increments of time are determined by the motion of 
the system via ephemeris time 

5t = N6X (6.8) 

where 



x = \?£r¥k (6-9) 



2{E-V) 



that has a corresponding expression in phase space 



N=J V * bpaP \ (6.10) 



2m(E - V) 

For the dynamical evolution of some / = f(q a Pa) we have 

6f(q aPa ) = NSX{f,H} (6.11) 

= NSX{f,^^ + U-E} (6.12) 

which is a physical motion! Now in analogy to what we previously did, we have 

Af = {f,H}(N5X) (6.13) 

And by changing N we inadvertently change A because ephemeris time determines the parametrization. So if we are 
to consider changing N we would change <5A, or vice versa, in such a manner that there is no diference and now there 
is no ambiguity in the evolution of /. We see this as a clear indication to choose a phase space with no preferred 
apriori assumed parametrization and the presymplectic manifold is just what we need. 



VII. PHASE SPACE AS A PRESYMPLECTIC MANIFOLD 

We recall from the symplectic formulation that the Legendre transformation is the fibre derivative that takes us 
from the configuration to the phase space. 

For some Lagrangian L(q, (q)) e C the fibre derivative is a map FL:TC — > T*C In the usual unconstrained case 

PoFL=™ (7.1) 

exists for all q's and so 



H o FL — pq - L 



(7.2) 
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is unconstrained. Thus we attain a symplectic form on phase space which is strongly non-degenerate, i.e there exists 
an X H e TT*C such that 

lx h u = dH (7-3) 

This is due to Lxw=0 that implies Lxdw+dixw=0 sinceT*C is simply connected. This gives Lxuj=dH and X=Xh- 
Thus, in such a system where the fibre derivative is a total diffeomorphism, the map 

b : TT*C -> T*T*C (7.4) 

is bijective. Thus the symplectic two-form is strongly non-degenerate. We shall refer to this as the symplectic case. 
Henceforth for simplicity we shall write T*C = T. In the completely constrained case given in our example above 

H = v a( f> a (7.5) 

the fibre derivative H o FL = giving a singular Lagrangian. 

Thus the condition that the symplectic form be strongly non-degenerate is to be relaxed and so we choose a weakly 
non-degenerate presymplectic form. The question that now arises is as to where and how exactly does the map b act as 
a bijection (locally)? The answer is that it acts on the foliations of phase space which are the constraint submanifolds 
and it does so by means of the presymplectic algorithm. 



VIII. PRESYMPLECTIC FORMALISM 

On the total phase space given by (T,uj) the symplectic case gives 

l x uj = dH (8.1) 

but here, in the first class case (it can easily be generalized into any number of constraints, but we choose first class 
for simplicity) there exists a constraint submanifold r C V on which the two form is strongly non-degenerate. Thus 
we choose an inclusion mapping 

h ■ f -> r (8.2) 

for which 

(t x uj - dH) o ji = (8.3) 
In general, the presymplectic algorithm generates a string of submanifolds, 

A- f 3 f 2 f -A r (8.4) 

and this ends up in three cases- 
Case 1: We hit an Tj such that T; = $. This plainly means that the manifold does not describe a dynamical system. 
Case 2: We get a T/ such that dim(rj) = 0, Ti ^ <I>. This means that the system is static, i.e X = 
Case 3: There exists an I such that T; = Ti+i. Thus we have a whole set of submanifolds with dim(Ti) ^ 0. This 
is the only case that gives us physically relevant information. And from the equivalence of this algorithm and the 
Dirac-Bergmann theory of constraints, this is the submanifold on which all the constraints are satisfied. Thus we can 
compound 

]n ° jn-1 ° ••• ° h ° 32 ° j = 11 ( 8 - 5 ) 

and so, 

n : f „ — > T (8.6) 
Putting it all together we find that for completely constrained systems 

[i n , x uj - dH] t(n) = (8.7) 

And the only condition T(n) needs to satisfy is 

T p T n = ker(cu),peT(n) (8-8) 
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This characterizes the arbitrariness in the foliation. We shall now se its implementation for completely constrained 
systems. We introduce some notation that we shall use throughout the paper. Let us, for simplicity, consider a 
completely constrained first-class Hamiltonian system. 

H = tffa = (8.9) 

where the Vi are arbitrary Lagrange multipliers and fa are the constraints. 
The equation 

gives us the constraint 

fa(q a ,p a )=Q (8.11) 

which defines T. We now choose a local coordinate patch and define 

uj = dp a Adq a (8.12) 

The kernel of this will be the Hamiltonian vector fields whose integral surfaces define the motion of the system entirely. 
Thus the aim of using the presymplectic equation is to obtain the Hamiltonian vector field. And so 

[l^xuj - dH\r M = (8.13) 

gives us tt*X = Xh- now, when we take into account the vanishing of the constraints on the sub manifold, we attain 
the actual geometric form of the presymplectic equation 

™ b |r (n) =0 (8.14) 

(henceforth just T will denote the final constraint sub manifold) We now define the flow of a function along these 
integral surfaces as the map 

fh[g(q,p)} = 9(q,p^) (8.15) 

where A is the parameter arising from the parametrization of the flow oiXjj which is again, completely arbitrary.The 
flow is determined by solving the Cauchy problem 

f°H = 9 

±f x H [g] = X H [g] (8.16) 

which is solved by 

00 in 
n=0 U - 

Now Hamilton's equatons are defined as the evolution of the flow of 

z a = [ q% I (8.18) 



Therefore, 



giving, 



And so the solution of Hamilton's equations is 



j- x &[9] = XH[z a ] (8.19) 



f H [z a }= [ q ) (8.20) 



/^ a ]=E^r. (8.2i) 

„=o n - 
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1. Relation to the Hamilton- J acobi theory 

Consider, on T, the symplcctic form 

Tr*Lo = dpAdq = iT*d6 = d(Tr*9) (8.22) 
If we choose a local coordinate patch (q,p) on T for which 

6 = pdq 

d6 = uj = dpAdq (8.23) 



Assuming that the phase space is simply connected, 



d(-K*6-6) = 

ir*6-0 = dS (8.24) 



Pi 

p l (8.25) 



This gives 

dS 
dQi 

dqt 

Assuming the Hamiltonian constraint has the form 

(9 ab PaPb + V(q a )) — E = H (8.26) 

£ = gives 

g ab PaPb + V(q a ) = E (8.27) 



and from the above 



which is the Hamilton- Jacobi equation. 



9ab 7$^ + V(qa) = E (8 ' 28) 



A. Relational particle dynamics 

In this section relational particle dynamics, where we basically try and eliminate all the redundancy present in 
many particle dynamical theories. We begin by considering the configuration space for an n-particle dynamical 
system moving in a 3 dimensional space: 

C = R 3 " (8.29) 

now this configuration space has redundency by virtue of translational, and Rotational invariance, which constitute 
the actions of the Euclidean group Euc(3) which will have 6 dimensions (3 rotations + 3 translations), so we shall 
consider the Relative Configuration space: 

C re i = R 3n /Euc(3) (8.30) 

Is this sufficient? - if we were considering a a universe consisting of but 3 particles, such that the only meaningful 
notion of position is the separation between the particles , we sould find that for the size of the triangles or straight 
lines formed by these particles the area and length should have no meaning whatsoever, this is due to the fact that a 
truly relational system isn't embedded in any background space, and so scale too must be eliminated, thus in totality 
we have to factor out scale as well, leaving us with a 3N-7 dimensional "shape space". The phase space will be the 
manifold T*C 3 (q a ,p a ) and 

q a = b(q a + c a - e abc d b q c ) 
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is the "corrected" position. The b,c and d are the auxialry gauge variables corresponding to dilations (scale changes), 
translations and rotations respectively. Now, the Hamiltonian for this theory is given by 

^ = ^(E|£ + fo2 ^))-£) ( 8 - 31 ) 

a 

+v 1 (J25 ab p a q b )+v 2 (J2Pa)+vUj2 eabC 1bPc) (8-32) 

Firstly, its evident that this theory is an example of a completely constrained Hamiltonian system having only first 
class constraints, secondly we find that the auxiliary gauge variables in the corrected position are actually the velocities 
with respect to the lapses of each of the parts of the Hamiltonian. This Hamiltonian is of the form 

H = NH + Y,Va<Pa 



Now, we can attain the Presymplectic equation by the following steps: 


first, we find that 




§-<>-> + 

a 


= E 


(8.33) 


^i = o=>D = y25 ab p aqb 

Sv a ^ 


= 0; 


(8.34) 


T a = J2p a 


-0; 


(8.35) 


ja = Y^e abc q bPc 


= 0; 


(8.36) 



The above constraints are known as the " Jacobi Constraints" . Now the Jacobi constraints define f and we find that 

f = T*R 3n /Euc(3) x V 

where R 3 ™ / Euc(3) x V is Shape space (D is the group of dilations) . We can move forward and attain the Presymplectic 
equation. 

[ix H (dp a Adq a )-dH]\ f =0 (8.37) 

which gives us: 

X n = {N— + v 2 +v iqa + (8.38) 
m a 

(8-39) 

-i Nb2 ^ - V ^ + £acb «3>c)^ (8.40) 
dq a dp a 

which satisfies: 

P^rif=0 (8.41) 
Now we can identify what part of the Hamiltonian constraint is responsible for physical evolution: 

£ ™ - lN ^i - (8 ' 42) 

(i.e. the vector field generated by the Hamiltonian constraint). It isn't hard to see that 

p a = m£ NH [q a ] 
and so, by plugging this back into |^ = we find that: 



\ m£ NH [q a ]£ NH [q a ] 
V 2(E-U) 1 ■ ' 



Thus we have re-derived the Ephcmcris-Lapsc correspondence. 
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B. The Relativistic Free Particle 

Now the above formalism will be applied to the relativistic free particle traveling through Riemannian (or Pseudo 
Riemannian) spacetime. Note that this chapter is presented purely to exhibit some of the features of presymplcctic 
dynamics, but not as a gravitational theory, for which we prescribe to the three-space approach (in the next paper). 
As of now, we shall consider the Hamiltonian 

H = N{g^ PlxPv - m 2 ) (8.44) 

and we arrive at T via: 

8T-L 

— = => g^v^Vv = m 2 

where g 1 *" is a Riemannian (or Lorentzian) metric. This defines a mass-hyperboloid fC m and we identify this with 
f 3 (p^,q^). As we did in the previous section, we attain the Hamiltonian vector field tautologically via 

[iXuidPnAdq^-dHWf = (8.45) 

and in this case we get 

Xn=N ^-- N J- p ^- (8 . 46 ) 



and it satisfies 



(X n ) b k m = (8.47) 



. By the very same calculation which we did for the previous system, we can find the Lapse, and here it turns out to 
be 



. where w M = X H [q^]. 



N = J 9 -^ (8.48) 



C. Geodesies and Hamiltonian Flows 



In the previous section, the Hamiltonian vector field satisfying the presymplectic equation for a free relativistic 
particle was obtained now we shall attain the solutions via the flow of this vector field. In particular let us condicr 
the flow of the position with respect to this Hamiltonian vector field, i.e 

n— oo 

= E -^ni^} (8-49) 



n! 

n=0 



from this, we know that 



now X n [q»] =pf*/N so, 



By the definition of a flow/^ [/] = g M (A), and putting all of this back into equation (8.50), and noting that 
g^v^Vv = g^v^v" we get 

<9V N c 

The above equation (for on shell particles) is the geodesic equation. This shows us that the Hamiltonian flows on 
phase space and the geodesies on configuration space physically describe the very same dynamics of the system. 
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IX. CONCLUDING REMARKS 

In this paper it has has been shown that the presymplectic equation is adequate to describe the dynamics of 
reparameterization invariant, completely constrained Hamiltonian systems. In all the systems considered thus far, 
the constraints saturate the theories, and even thought they help define the specific domain of their dynamics on 
phase space i.e. the constraint submanifold via the presymplectic algorithm, they and so the total Hamiltonian, do 
not explicitly feature in the final equation. The Hamiltonian vector field is only derived from the natural tautology. 
Also, the novelty of our formalism lies in the adoption of the interpretation of Barbour et. al. i.e. Machian Temporal 
Relationalism and its incorporation into a canonical phase space- framework, which can also account for the Ephemeris 
time and the violation of Dirac's theorem with regard to first class Hamiltonian constraints. More shall be said about 
the Barbour interpretation in our next paper in the context of canonical General Relativity. 
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